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ABSTRACT 

In this paper we report on an XMM-Newton observation of the ultraluminous infrared 
QSO Mrk 1014. The X-ray observation reveals a power-law dominated (T w 2.2) spec- 
trum with a slight excess in the soft energy range. AGN and starburst emission models 
fit the soft excess emission equally well, however, the most plausible explanation is 
an AGN component as the starburst model parameter, temperature and luminosity, 
appear physically unrealistic. The mean luminosity of Mrk 1014 is about 2 x 10 
erg s _1 . We have also observed excess emission at energies greater than 5 keV. This 
feature could be attributed to a broadened and redshifted iron complex, but deeper 
observations are required to constrain its origin. The light curve shows small scale 
variability over the ~11 ks observation. There is no evidence of intrinsic absorption in 
Mrk 1014. The X-ray observations support the notion of an AGN dominated central 
engine. We establish the need for a longer observation to constrain more precisely the 
nature of the X-ray components. 
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1 INTRODUCTION 



Mrk 1014 (V = 15.7; z = 0.163), with a far-infrared (FIR) 
, luminosity in excess of 1O 12 L (Yun 2001), is one of the 
5_j ■ brightest members of the class of 'warm' ultraluminous in- 
' frared galaxies (Sanders et al. 1988a) which from their warm 
IRAS 25/im/60/im colors and optical spectra, are believed 
to host powerful AGN. Sanders et al. (1988b) report on the 
detection of CO (1-0) emission and suggest that Mrk 1014 
may be an important link in the evolution of ultraluminous 
infrared galaxies into UV-excess quasars. 

Spectroscopy using the Infrared Space Observatory 
(ISO) adds quantitative information to this picture. The low 
resolution ISOPHOT-S spectrum presented by Rigopoulou 
et al. (1999) is dominated by AGN-like continuum emission, 
with no evidence for the aromatic 'PAH' emission typical 
for starbursts. The ratio of 7.7/im PAH feature and local 
continuum is <0.63, well below starburst values of ~3. The 
diagnostics proposed by Genzel et al. (1998) and Tran et al. 
(2001) then suggest an upper limit of «25% for the star- 
burst contribution to the infrared luminosity of Mrk 1014. 
Sturm et al. (2002) present ISO-SWS mid-IR fine structure 
line data for Mrk 1014. The source is at the limit of the the 
SWS sensitivity, and only tentative detections are obtained 
for the low excitation [NeII]12.8^m line and the high exci- 



tation [OIV]25.9/im line, with a ratio consistent with AGN 
dominance. 

Hubble Space Telescope (HST) NICMOS observations 
(Scoville et al. 2000) reveal twisting spiral isophotes beneath 
the dominant QSO nucleus, indicating either a starburst spi- 
ral disk or tidal debris. Optical observations of the host 
galaxy portray a prominent tidal arm with several knots 
(Surace & Sanders 2000; Surace et al. 1998). The tidal arm 
is composed primarily of intermediate age stars (~ 1 Gyr) 
with very little contribution from older stars (Canalizo & 
Stockton 2000). Nolan et al. (2001) determine that the ma- 
jority of the flux is associated with an instantaneous burst of 
star formation at approximately 12 Gyr, however, Canalizo 
& Stockton (2000) find several regions of recent (~ 0.2 Gyr) 
star formation accounting for up to 30% of the total lumi- 
nous mass along the line of sight. Both age estimations are 
too high to be associated with on-going starburst activity 
(~10 Myr), which produces X-ray emitting hot gas. 

Mrk 1014 was observed with XMM-Newton (Jansen et 
al. 2001) as part of the guaranteed time programme. The ob- 
jective of the study was to obtain a high quality X-ray spec- 
trum of an AGN dominated ultraluminous infrared galaxy, 
and determine if starburst activity is a necessary component 
to explain the spectral energy distribution. 

In the following section we will present the observa- 
tions and discuss the data reduction. In section 3, we will 
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describe the spectral models and discuss their significance. 
A description of the light curve will be presented in section 
4. In section 5 we will summarize our findings. A value of 
the Hubble constant of Ho=70 km s _1 Mpc -1 and a cosmo- 
logical deceleration parameter of qo = 5 have been adopted 
throughout. 



2 OBSERVATION AND DATA REDUCTION 

Mrk 1014 was observed with XMM-Newton on 2000 July 
29 during revolution 0117 for a duration of 15 ks. All 
instruments were functioning normally during this time. 
The European Photon Imaging Camera (EPIC) pn detec- 
tor (Striider et al. 2001) was operated in full-frame mode, 
and the two MOS cameras (Turner et al. 2001) were oper- 
ated in large-window mode. All three EPIC CCD cameras 
used the medium filter. 

The MOS and pn Observation Data Files (ODFs) were 
processed to produce calibrated event lists using the XMM- 
Newton Science Analysis System (SAS) v5.3. Unwanted hot. 
dead or flickering pixels were removed as were events due to 
electronic noise. Event energies were corrected for charge- 
transfer losses. The latest available calibration files^| were 
used in the processing. Single and double events were se- 
lected. Background light curves were extracted from these 
event lists to search for periods of flaring, but showed the 
background to be stable throughout the duration of the ob- 
servations. The total amount of "good" exposure time was 
10534 s with the MOS cameras and 11298 s with the pn 
camera. 

The source plus background photons were extracted 
from a circular region of radius 40 arcsec for the pn and 
MOS data. The centroid position of Mrk 1014 in the EPIC 
images is a=0l' l 59 m 50 s .2, <5=+00°23'41.2", which is in good 
agreement with the optical position given by Argyle & El- 
dridge (1990). Background events were extracted from an 
off-source circular region roughly 90 arcsec in radius. 

The Galactic column density toward Mrk 1014 is (2.6 
± 0.2) x 10 20 cm" 2 (Dickey & Lockman 1990). 



3 SPECTRAL ANALYSIS 

The source spectra were grouped such that each spectral bin 
contains at least 60 counts for the pn data and 20 counts 
for the MOS. Spectral fitting was accomplished with the 
XSPEC vll.l software package (Arnaud 1996). The quoted 
errors on the derived best-fitting model parameters corre- 
spond to a 90% confidence level unless stated otherwise. 
When we calculate luminosities, we shall assume isotropic 
emission unless stated otherwise. 

Displayed in Figure 1 is the EPIC pn spectrum along 
with the background spectrum in the 0.3-10 keV energy 
range. The background is no longer negligible at energies 
higher than about 8 keV. The same is true for the MOS 
observations. We will only consider the 0.3-8.0 keV spectral 
range when generating models of Mrk 1014. In the initial 

* ml_mcdv9q20t5r6_all_15.rsp and 
m2_medv9q20t5r6.all_15.rsp for MOS1 and MOS2, respectively 
and epn jf20_sdY9_medium.rmf for the pn. 
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Figure 1. The dominant spectrum is the source+background ob- 
served with the EPIC pn camera. The weaker spectrum, averaging 
approximately 0.01 counts s _1 keV -1 , is the background contri- 
bution. The background spectrum becomes just as significant as 
the source spectrum at energies higher than 8 keV. The same is 
true for the MOS observations. 

models the photoelectric absorption was treated as a free pa- 
rameter. Absorption above the Galactic Nh was not found 
in any spectral model indicating very little intrinsic absorp- 
tion in Mrk 1014. In the remainder of the analysis we have 
fixed the absorption parameter to the Galactic value stated 
above. 

A simple power-law modified by Galactic absorption 
was fitted separately to the pn and combined MOS 0.3-8 
keV spectra with mediocre results (x 2 = 160.4/159 dof for 
the pn data, and % 2 = 252.6/229 dof for the combined MOS 
data). As we will demonstrate in section 3.1, the addition of 
a soft X-ray component significantly improves the fits. The 
photon indices were V = 2.46 ± 0.03 and F = 2.42 ± 0.04 for 
the pn and MOS data, respectively. Despite the reasonable 
agreement between the pn and MOS data we do not combine 
the spectra in order that we may search for discrepancies 
in the more complicated models discussed below. We will 
alternate our discussion between the two spectra and will 
explicitly state any incongruity between models. 

There is no convincing evidence of a narrow iron Ka 
emission in either spectra nor is there evidence of spectral 
variability. A hardness ratio light curve was produced using 
the 0.3-2 keV and 2-8 keV energy bands to test for spectral 
variability, but the hardness ratio light curve is represented 
by a constant fit (\ 2 ~ 23.5/55 dof). 

3.1 The Soft Energy Spectrum 

The statistics indicate that the fits are significantly im- 
proved by the addition of a soft component to the spectral 
models. Given the nature of Mrk 1014, this soft component 
may be attributed either to AGN or possible starburst ac- 
tivity in the nucleus. We put forth a number of models to 
test both hypotheses. A simple blackbody, double blackbody 
and multiple blackbody disc models are used to search for 
an AGN contribution; thermal bremsstrahlung and a Ray- 
mond & Smith (1977) models are suggested to investigate 
the starburst nature. In the remainder of this section we 
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Table 1. Spectral Parameters of the Attempted EPIC pn Models 



Figure 2. The 0.3-8 keV EPIC pn and combined MOS spectra 
of Mrk 1014 and model components are displayed in the upper 
panel. The data are fitted with a blackbody plus power-law com- 
ponents (x 2 = 145.7/156 dof and 210.7/212 dof for for the pn 
and MOS, respectively). The solid line is the pn model and the 
broken line is the combined MOS fit. The lower panel show the 
residuals (data/model) of the models. Notice the excess residuals 
at energies higher than 5 keV. 



will discuss these models. We have summarized the results 
in Table 1. 

As can be seen in Figure 2, a blackbody plus power-law 
provides an acceptable fit to the data (x 2 = 145.7/156 dof 
for the pn data, and \ 2 ~ 210.7/212 dof for the combined 
MOS data). The derived power-law indices are similar in the 
pn model (F = 2.28 ± 0.06) and the MOS model (r = 2.21 
± 0.08). The unabsorbed luminosities for the model com- 
ponents are also consistent: We measure a total 0.3-8 keV 
luminosity of 1.65 x 10 44 erg s _1 with MOS, of which the 
blackbody component accounts for 16%. With the pn model 



we calculate a luminosity of 1.58 x 10 44 erg s" 1 , of which 
13% is attributed to the blackbody component. A double 
blackbody plus power-law also provided a good fit for the 
pn (x 2 = 139.5/154 dof) and MOS {x 2 = 209.9/210 dof) 
data, however several parameters were unconstrained, such 
as the temperature of the second blackbody. This would in- 
dicate that the model is not an accurate representation of 
the data given the available photon statistics. A multiple 
blackbody disc in addition to the power-law was also a rea- 
sonable fit to the data (x 2 = 142.2/156 dof and 210.3/212 
dof for the pn and MOS, respectively). 

In an attempt to model a possible starburst component 
we employed the Raymond-Smith model for an emission 
spectrum from a hot, diffuse gas. In combination with a sim- 
ple power-law the models fit the data well (x 2 = 152.2/156 
dof for the pn data and x 2 = 215.8/212 dof for the MOS 
data). We find kT = 273±5°° eV, which is lower than the 
temperature of >600 eV usually found in X-ray starbursts 
(Ptak et al. 1999). This claim is further supported by the 
recent XMM-Newton observations of the ULIRG NGC 6240. 
The starburst in NGC 6240 is adequately fitted by two 
Raymond-Smith components with temperatures kT = (0.62 
± 0.20) keV and (1.11 ± 0.10) keV (Keil, Boiler & Fuji- 
moto, in preparation). We calculate an X-ray luminosity of 
~ 5 x 10 42 erg s _1 for the Raymond-Smith component in 
Mrk 1014. This X-ray luminosity appears to be surprisingly 
large for a starburst, therefore, the detection of a starburst 



Soft Energy Model 


kT ( 


eV) 


Power-law Index (r) 




Blackbody 


147 ± 15 


2.3 ± 0.1 


0.93 


Double Blackbody 1 


74 +27 


n i+0.2 


0.91 


Mult. Blackbody Disc 2 


179" 


his 

-9 


2.2 ± 0.1 


0.91 


Raymond-Smith 


273 4 


-100 
-50 


2.4 ± 0.4 


0.98 



x Data is presented for only one of the blackbodies as the second 
blackbody was not constrained. 

2 The temperature presented is that at the inner disc radius. 



Table 2. The 0.3-8 keV Unabsorbed Luminosity of the Various 
Model Components (xlO 44 erg s — 1 ) 



Soft Energy Model 


Lsoft 


^power — law 


Lsoft 
L total 


Blackbody 


0.20 


1.38 


0.13 


Double Blackbody 


0.45 


1.16 


0.28 


Mult. Blackbody Disc 


0.30 


1.29 


0.19 


Raymond-Smith 


0.05 


1.55 


0.03 



component above the luminous AGN power-law is very un- 
likely. A second respectable fit was obtained by adding a 
thermal bremsstrahlung component to the power-law (x 2 = 
139.1/156 dof, 210.1/212 dof for the pn and MOS, respec- 
tively) , however the resulting temperatures and luminosities 
are unphysically large. 

On a statistical bases, none of the attempted fits appear 
superior over the others. However, from physical arguments 
the AGN dominated spectral model appear the most realis- 
tic explanation for the X-ray features of Mrk 1014. In Table 
1 we present the parameters of the various models we at- 
tempted for the pn spectrum. Results are comparable for 
the MOS observations. In Table 2 we have calculated the 
unabsorbed luminosities measured for the soft and power- 
law component for each of the pn models used. Again, the 
results are comparable with the MOS observations. 

The spectral analysis reveals that Mrk 1014 is AGN 
dominated. This is further supported by comparing the soft 
X-ray (0.1-2.4 keV) flux to the FIR (40-120 fim) flux. Boiler 
& Bertoldi (1996) found that in equilibrium, the ratio of the 
X-ray to FIR flux is about 10~ 2,5 . A flux ratio less than 
10 -2 ' 5 indicates that starburst activity is dominant. The 
total FIR flux of Mrk 1014 is 9.93 x 10" 11 erg s _1 cm -2 , 
computed from the IRAS 60 and 100 /im fluxes following 
Helou et al. (1985). To compute the soft X-ray fluxes we used 
the models described in section 3. The extrapolated 0.1-2.4 
keV XMM-Newton flux is 3.68 x 10" 12 erg s" 1 cm" 2 . We 
calculated a value of 10 -1 ' 4 for the flux ratio in Mrk 1014, 
an order of magnitude greater than what is expected from 
a starburst. 



3.2 The Hard Energy Excess 

As can be seen in Figure 2, there appears to be a gentle 
uprising in the residuals at energies higher than 5 keV. This 
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reveals a power-law dominated spectrum which may be up to 
thirty times more luminous than the soft X-ray excess. The 

i 1 1 1 1 1 1 1 1 1 1 1 1 existence of an X-ray starburst is highly unlikely due to the 

low temperature and high luminosity found in the models. 
Other evidence, such as a variable light curve and the soft X- 
ray to FIR flux ratio suggest an AGN powered continuum in 
Mrk 1014. A longer observation with XMM-Newton would 
allow us to better constrain the nature of the soft energy 
component. In addition, a longer look will allow us to probe 
the interesting high energy excess above 5 keV and deter- 
mine if it is connected to a broadened and redshifted iron 
complex. Clearly, more multi-frequency observations are re- 
quired to disentangle and compare the starburst and AGN 

1,1,1,1,1,1, contributions in different energy bands. 

2000 4000 6000 8000 10000 12000 

Time (s) 



Figure 3. The EPIC pn light curve of Mrk 1014 in the 0.3-8 keV 
energy range with a bin size of 400 s. The ordinate is the count 
rate displayed in count s — 1 , and the abscissa is the time offset 
from the start of the exposure on 2000 July 29 at 21:29 UT. 

trend is apparent in all of the models presented. By adding a 
second power-law with a negative slope to the models we find 
an improvement in the fits. Using the F-test for the addition 
of one free parameter, we find that we have a significant im- 
provement (>99.95 %) in the quality of the fits. The current 
observation does not allow us to constrain the nature of this 
excess high energy emission. It would be especially interest- 
ing to determine if this hard excess can be attributed to a 
broadened and redshifted iron complex. Longer observations 
are required to further investigate the characteristics of this 
hard energy excess. 



4 TIMING ANALYSIS 

In Figure 3 we present the 0.3-8 keV EPIC pn light curve in 
400 s bins. The MOS light curves are similar to the pn light 
curve, but only the pn data will be presented here given the 
higher signal to noise. 

The light curve of Mrk 1014 does not show much in the 
way of interesting features but there does appear to be some 
small scale variability during the short 11 ks observation. 
A constant model through the mean of the light curve is 
not an acceptable fit (x 2 = 39.1/27 dof). The mean pn 
count rate during the observation is (1.18 ± 0.09) counts 
s _1 . There is a difference of less than 0.3 counts s _1 between 
the strongest and weakest events in the pn light curve. The 
largest fluctuations in the light curve are observed at ~1600 
s and ~8000 s, in which the count rate increases by 0.2 
counts s _1 in 400 s. The limited statistics in these peaks do 
not allow us to ascertain more information. The mean count 
rate measured with the MOS data is 1.17 ± 0.07 counts s _1 . 



5 SUMMARY 

The 0.3-8 keV spectrum of the Mrk 1014 was fitted with a 
number of models to determine if the soft energy emission 
above a predominant power-law is driven by an AGN central 
engine or starburst activity. The XMM-Newton observation 
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